Confronting Lemaitre-Tolman-Bondi models with 
Observational Cosmology 



Juan Garcia-Bellido 1,2 , Troels Haugb0lle 13 

1 Instituto de Ffsica Teorica UAM-CSIC, Universidad Autonoma de Madrid, 
Cantoblanco, 28049 Madrid, Spain, 

2 Kavli Institute for Theoretical Physics, University of California Santa Barbara, CA 
93106-4030, 

3 Department of Physics and Astronomy, University of Aarhus, DK-8000 Aarhus C, 
Denmark 

E-mail: juan.garciabellido@uam.es, haugboel@phys.au.dk 

Abstract. The possibility that we live in a special place in the universe, close to the 
centre of a large void, seems an appealing alternative to the prevailing interpretation 
of the acceleration of the universe in terms of a ACDM model with a dominant dark 
energy component. In this paper we confront the asymptotically flat Lemaitre-Tolman- 
Bondi (LTB) models with a series of observations, from Type la Supernovae to Cosmic 
Microwave Background and Baryon Acoustic Oscillations data. We propose two con- 
crete LTB models describing a local void in which the only arbitrary functions are the 
radial dependence of the matter density SIm and the Hubble expansion rate H . We find 
that all observations can be accommodated within 1 sigma, for our models with 4 or 5 
independent parameters. The best fit models have a \ 2 very close to that of the ACDM 
model. A general Fortran program for comparing LTB models with cosmological ob- 
servations, that has been used to make the parameter scan in this paper, is made pub- 
lic, and can be downloaded at http : //www . phys . au . dk/~haugboel/sof tware . shtml 
together with IDL routines for creating the Likelihood plots. We perform a simple 
Bayesian analysis and show that one cannot exclude the hypothesis that we live within 
a large local void of an otherwise Einstein-de Sitter model. 
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1. Introduction 

Cosmology has been traditionally prone to speculation, a science which has had more of 
wishful thinking than actual deduction. Even today, with the extraordinary revolution 
in observational cosmology, there are many assumptions we have to take for granted, in 
order to interpret the present observations. In particular, the present abundance of less 
that 5% of critical density in matter we know and measure in the laboratory, i.e. the 
baryons we are made of, compared to the 23% of unknown "dark matter" and the 72% 
of even more unknown "dark energy", seems a rather peculiar and perhaps suspicious 
composition for the universe, given their disparate evolution rates. Are we sure we have 
now the correct picture of the universe? 

We live in a very isolated part of the universe and have not reached that much 
beyond our solar system. All we can say about the universe has been inferred from 
observations done under very general assumptions; but some of them, no matter how 
fundamental they may seem, may be wrong. It is worth pointing out how Einstein failed 
to make one of the most fundamental predictions of his new theory of general relativity, 
and had to introduce an "absolute" in the theory — the cosmological constant — much 
to his regret, because he had the prejudice, supported by incomplete observational data, 
that we lived in a static universe without a beginning or an end. Soon after the discovery 
of the "recession" of galaxies by Hubble did he renounce to his prejudice on the validity 
of the Perfect Cosmological Principle, which assumed maximally symmetric space- 
times. For 70 years cosmologists have worked under the assumption of the less strict 
Cosmological Principle, which imposes maximal symmetry (homogeneity and isotropy) 
only for the spatial sections. But are these symmetries consistent with observations? It 
is evident to anyone that looks at the sky in a clear night far away from city lights that 
the universe is not homogeneous and isotropic. 

It has usually been argued that these fundamental symmetries should only be 
"expected to be valid on very large scales", but how large are large scales? In fact, 
the distribution of matter in our local vicinity, i.e. within several Mpc, is very far from 
homogeneous; density contrasts reach enormous values not only at the centres of galaxies 
but also on larger scales like clusters and superclusters, stretching over hundreds of Mpc. 
So, how far do we have to go before we reach homogeneity? Certainly present galaxy 
catalogues are not really there yet, and it has been speculated that perhaps with the 
next generation of deep catalogues like DES [1] we will finally reach the homogeneity 
limit. But, if we do not live in a homogeneous universe, how do we interpret the 
observations we have of objects whose light has travelled a significant fraction of the 
age of the universe in order to reach us? Our present assumption is that, in practice, 
the intervening inhomogeneity averages out and everything works as if we lived in a 
homogeneous universe. For several decades this assumption has been a valid one, and has 
provided confidence into the construction of the so-called Standard Model of Cosmology. 

It is only now that we begin to have sufficiently good cosmological data, and 
certainly will have even better data in the near future, that we can be critical 
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and pose the appropriate questions. The issue of homogeneity of the universe has 
been often dismissed because of the apparent extraordinary isotropy of the cosmic 
microwave background. However, any mathematician can readily show that isotropy 
and homogeneity are very different sets of symmetries and one does not imply the other. 
Nevertheless, if we impose the further assumption, usually stated as the Copernican 
Principle, that any point in space should be equivalent to any other, i.e. that we don't 
live in any special place in the universe, then a mathematical theorem states that if all 
equal observers see the universe isotropic around them, then the universe must be not 
only isotropic but also homogeneous. What remains to be proven is that all observers 
are in fact equivalent in the patch of the universe we call the observable universe. Some 
will be near a large concentration of mass and others will be in large voids. Certainly 
what these two types of observers see will differ from what an idealised observer living 
in a perfectly homogeneous universe would see. Unfortunately, we have never spoken to 
anyone at the other side of the universe. 

The advantage of the present state of Cosmology is that we can begin to pose those 
questions and hope to get concrete answers, while just ten years ago it would have been 
futile. Moreover, like with the first inclusion of the cosmological constant in the theory, 
almost 100 years ago, the physics community is very much puzzled about the nature of 
this so-called vacuum energy. Its properties defy our basic understanding of quantum 
physics and, moreover, it reminds us suspiciously of the Maxwellian ether, which led 
the way (via its disappearance) to a new understanding of physical reality. To put 
the question straight, are we sure we live in an accelerated universe which is driven by 
some unknown vacuum energy? Could it be that we have misinterpreted our superb 
cosmological data and what those photons coming from afar are telling us is something 
completely different? 

The last few years we have seen a tremendous burst of activity, both at the 
theoretical and observational level, in order to disentangle the subtle intricacies of the 
actual data sets from their interpretation. As with any hard science (and Cosmology 
is indeed finally becoming one), such enterprises can be approached only via further 
observational crosschecks. It is no longer true that "astrophysicists are often in error, 
but never in doubt", as Lev Landau once said to annoy his colleague Yakov Zel'dovich 
[2]. We can now propose new ways to measure more observables in a wider theoretical 
construction. Perhaps it is time to explore the troubled waters of non-maximally 
symmetric spatial sections of the universe. In particular, since we indeed observe a 
high degree of isotropy in the cosmic microwave background, we can start by exploring 
the simplified version of a spherically symmetric inhomogeneous model, which comes 
under the name of the Lemaitre-Tolman-Bondi (LTB) model [3j HJ [5]. 

While ordinary Friedmann- Robertson- Walker (FRW) space-times are characterised 
by two functions, the Hubble rate H(t) and the density parameter Q(t), which depend 
on cosmic time but are independent of the radial coordinate, the LTB models have also 
two arbitrary functions, H(r,t) and Q(r, t), which depend on both time and the radial 
coordinate. The Einstein equations for the LTB model are sufficiently simple that they 
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can be integrated exactly in terms of two arbitrary boundary conditions. Thus, the LTB 
models have lower symmetries, and thus more freedom, but nevertheless make specific 
predictions about the behaviour of light along geodesies of the new metric. Therefore, 
one can evaluate the corresponding observables, like cosmological distances (or ratios 
of distances), from the available sets of data. Surprisingly, as we will show in the next 
sections, present cosmological data does not yet seem able to exclude with confidence 
a universe which is not exactly homogeneous. This should not come as a surprise, we 
do live in an (locally) inhomogeneous universe, some say that within a large underdense 
void [HI El U M, EU similar to that which induces a cold spot in the CMB [T21 EH E] 
so why should we assume global homogeneity? In fact, the possibility that we happen 
to live in the centre of the world was advocated long ago in Ref. [TS], based on the 
stochastic inflation formalism [T6l [T7] . What is surprising is that within this framework 
one can account for (almost) all observational evidence without having to introduce 
an unknown in the theory, an absolute, whose properties are highly mysterious, not to 
mention the alternative highly artificial and ad hoc modifications of gravity on the very 
large scales. 

In this paper we propose an inhomogeneous model of the universe, with a local 
void size of a few Gpc and asymptotically Einstein-de Sitter, and then use the present 
sets of data (CMB, LSS, BAO, SNIa, HST, Ages, etc.) to constrain its parameters. In 
Section 2 we describe the general LTB metric, the Einstein equations and the definitions 
of cosmic distances. We also give a novel series expansion that allows one to integrate 
to arbitrary precision the Einstein equations for arbitrary functions Ho(r) and Qm^)- 
We then describe a specific model with a concrete form of these two functions that are 
plausible but simple matches to observations, and show the different (longitudinal and 
transverse) rates of expansion and the apparent acceleration. We also show that the 
comoving size of the sound horizon depends on the distance to the centre of the void 
in the LTB model, and relate it to the expansion rates at different redshifts. This is a 
prerequisite for making a realistic comparison with Baryon Acoustic Oscillation data. 
In Sect. 3 we present the data sets we use to constrain the model, with the subtleties 
needed to correctly interpret these in terms of the model. In Sect. 4 we give the main 
results and present a Bayesian analysis to ascertain the goodness of fit of the model by 
comparing it with a standard ACDM model with variable (but constant) equation of 
state parameter w. In Sect. 5 we give the conclusions. 

2. The LTB model 

The most general (cosmological) metric satisfying spherically symmetric spatial sections 
can be written as 



where dfl 2 = d6 2 + sin 2 9d(f> 2 . Assuming a spherically symmetric matter source, 




(1) 



= -p M (r,t)« 
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the (0, r) component of the Einstein equations, G° r = 0, implies X(r,t) = 
A'(r, t)/y/l — k(r), with an arbitrary function k(r) playing the role of the spatial 
curvature parameter. Note that we recover the FRW metric imposing the extra 
homogeneity conditions, 

A(r, t) = a(t) r , k(r) = kr 2 . 

The other components of Einstein equations read[18, 19J 
A 2 + k AA' k'ir) 

A 2 + 2AA + k(r) = . (3) 
Integrating the last equation, we get 

p f(0 _ Mr) m 
A 2 A 3 A 2 ' 1 ' 

with another arbitrary function F(r), playing the role of effective matter content, which 
substituted into the first equation gives 

Combining these two equations we arrive at 
2A lA' AttG 

which determines the "effective" acceleration in these inhomogeneous cosmologies. Note 
that the notion of acceleration becomes ambiguous since backward proper time has both 
a time and a spatial component. 

The boundary condition functions F(r) and k(r) are specified by the nature of the 
inhomogeneities through the local Hubble rate, the local matter density and the local 
spatial curvature, 

*™ = l$y (7) 

F(r) = H 2 (r)n M (r)A 3 (r) } (8) 
k(r) = H 2 (r)(n M (r)-l)A 2 (r), (9) 



where functions with subscripts correspond to present values, A (r) = A(r,to) and 
H (r) = H(r,t ). With these definitions, the r-dependent Hubble rate is written 
as I 



H 2 (r,t) = H 2 (r) 



An(r)\ 3 „ , „ f Mr" x 2 



MmTrH +(l-n M (r)) 



(10) 



A(r,t) ) y my " \A(r,t), 

Inhomogeneities come in two different classes: in the matter distribution or in the 
expansion rate, which are mutually independent. Moreover, the extra gauge freedom of 
the synchronous comoving gauge allows us to choose 

A(r,t ) = Ao(r) =r. (11) 
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Then we can integrate the Hamiltonian constraint equation (TTOT) to provide comoving 
time as a function of the radial coordinate, 

rA(r,t)/Ao(r) dx 

H {r)t{r) 



y/n M (r)/x + l-n M (r) 



A(r, t) j x + n M (r)A (r) _ Q M (r) Q K (r)A(r,t) ^ ^ 



A (r)y^(rj\l n K (r)A(r,t) ^{r) \j n M (r)A (r) 

where fi^r) = 1 — S^mO")- In particular, setting A(r,t) = A (r) we find the current 
age of the universe 

Ho(r)tHB(r) = 1 X^W. M= smh-sM. (13) 

For an observer located at the centre r = 0, by symmetry, incoming light travels 
along radial null geodesies, ds 2 = dQ 2 = 0, and time decreases when going away, 
dt/dr < 0, and we have 

dt_ = A'(r,t) 
dr y/l - k(r) 
which, together with the redshift equation, 
d\og{l + z) = A'(r,t) 
dr y/i - k(r) 

can be written as a parametric set of differential equations, with N = log(l + z) being 
the effective number of e-folds before the present time, 
dt A'(r,t) 



dN A'(r,t) 



(16) 



dr = y 7 ! ~ k(r) 
dN~ A'(r,t) ' 

from which the functions t(z) and r(z) can be obtained. From there one can immediately 
obtain both the luminosity distance, the comoving distance and the angular diameter 
distance function of redshift, 

d L (z) = (l + z) 2 A[r(z),t(z)], (18) 

d c (z) = {l + z)A[r(z) 1 t(z)] 1 (19) 

d A (z) = A[r(z),t(z)}, (20) 

2.1. Series solution 

In order to integrate out the redshift dependence it will be useful to make a series 
expansion of the cosmic time variable t(r) as a function of the space-dependent scale 
factor A(r,t). For this purpose we will define new variables 

V = H o(r)V^F) t(r) = \{5af' 2 , (21) 



Confronting Lemaitre-Tolman-Bondi models with Observational Cosmology 



7 



where S(r) is a generalised density contrast ratio and a(t) is the Einstein-de Sitter (FRW) 
scale factor, 



5(r)-. 
a(t) 



n M (r) 



2/3 



(23) 
(24) 



With these definitions, the time integral ffT2"]) can be written as y — v/ x(l + x) —\d\_\Jx- 
V 1 + x] , which can be expanded in series and inverted. With the definitions 



9W 

f(r) 

h(r) 



z + -z 
5 



:Z 6 + 



23 



1894 



175 7875 3031875 
H' {r) Q' M (r) fl + n M (r)/2 



z 5 + 0{z & ) , 



H (r) n M (r) V 1-JMr) 



1 

- + 



(25) 
(26) 
(27) 



r n M (r)(l-n M (r)) 
we can write the solution as a power series, whose coefficients can be calculated with 
arbitrary precision, 

(28) 

2 r 

3t5 
r 



A(r,t) = ^g(ad) 



A(r,t) - 
A'(r,t) 

A'(r,t) 

A(r,t) - 

A'(r,t) 



a 5 g'(a S) , 



8 



g (a5)h(r) + -a5g'(a5)f(r) 



2 r 
M5 



a5g'(a6)h(r) + ~a6g'(a5) f(r) + ^(a6) 2 g"(a5) f(r) 



F(r) 



2A 2 {r,t) ' 
F(r)A'(r,t) 



F'(r) 



(29) 
(30) 

(31) 

(32) 

(33) 



A 3 (r,t) 2A 2 (r,t)' 
These functions allow us to construct any other observable. For instance, the 
transverse and longitudinal rates of expansion can be written as 

A(r,t) 



H T (r,t) 
H L (r,t) = 



A{r,t) ' 
A'(r,t) 



(34) 



(35) 



A'(r, t) ' 

Note that in general these two functions will be different, and they enter into other 
observables. We can also construct quantities like the "effective" acceleration parameter 

d In H(z) 



q(z) 



•1 + 



d ln(l + z) 
where H(z) is in fact H L (r(z),t(z)). 



(36) 
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One could also define an effective equation of state parameter 



. . p(z) , 1 d ln 

W(z) = =-7-r = -1 + 



tf 2 (^)/tf 2 (r)-ft M (r)(l + z) 3 



(37) 



p(z) 3 dln(l + z) 

where H(z) is here Ht(t(z), t(z)). 

2.2. Parametric solution 

At fixed r the r-dependent Hubble rate Eq. (TlOT) is just like the normal Friedmann 
equation, and the standard way to explicitly solve for A(r, t) is to use an additional 
parameter rj. With the selected gauge the solution is 

^■«) = 2p^oll CMh M- 1 ^M (38) 

H ^ t= w^w lsiam - V] (39) 

Given r and t, by solving Eq. ( 1391) 77 can be found, and combining Eqs. ( l38i) -( [39j) with 
the Einstein equations we can then derive any necessary quantity. 



We have used both the series solution, implemented in a Mathematica notebook, 
and the parametric solution, implemented as a Fortran program, to make and double 
check all numerical computations in this paper. The Fortran 90 program together with 
a set of IDL routines for making likelihood plots is made publicly available and can be 
downloaded at htt p : //www . phys . au . dk/~haugboel/sof tware . shtml| 



2.3. The GBH model 

Here we define a new type of LTB model, which is completely specified by the matter 



content VLm{t) and the rate of expansion H (r), 

^)-«U + (^-«U)( 1 - t ^aS H )- <40) 

Ho{ ^H M+ ^-H^( i^_^ y (41) 

which is governed by 6 parameters, 

Q ont determined by asymptotic flatness (42) 

fij n determined by LSS observations (43) 

H out determined by CMB observations (44) 

H in determined by HST observations (45) 

r characterises the size of the void (46) 

Ar characterises the transition to uniformity (47) 
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d A (t BB ) [Gpc] d A (t BB ) [Gpc] 

Figure 1. The radial dependence of the physical matter density in units of the critical 
density (p m ), of our density function (f2 m ), that coincides with the usual matter density 
at the centre and asymptotically, and of the transverse and longitudinal expansion 
rates (Ht and Hi,). The radial axis is the angular diameter distance (dy4(r, £bb))j and 
everything is taken at the current time for the central observer. To the left (right) is 
shown the best fit GBH (constrained) model (see table [5]). Also shown to the right is 
an example of the density profile of a model with a sharp transition (Ar/ro = 0.3), 
but still within l-er of the best fit. 



We fix f2 out = 1 and let the other five parameters vary freely in our parameter scans 
(see tabled] for the priors). For instance, a plot of f2jv/(r), H T {r) and H L (r) for the two 
best fit models, as a function of the angular diameter distance today, d,A(r, £bb), can be 
seen in Fig. J. Also shown is the density profile of a model with a sharper transition 
(Ar/r = 0.3), but still within l-er of the best fit. This illustrates that observations 
allow for shallower density profiles close to the origin, and that for the physical matter 
density pm, we naturally get a shell-like transition. 



2-4- The constrained GBH model 



We have also considered a more constrained model, in which the Big Bang is 
homogeneous, that is, the spatial hypersurface at the Big Bang does not depend on 
the radial coordinate r. This can be obtained simply using Eq. ([TBI , by a choice of 
H (r), 



H (r) = H 



£l K (r) 



sinh 



( 4 * 



so that £bb = cHq 1 is universal, for all observers, irrespective of their spatial location. 
Note that in this model we have less freedom than in the previous model, since now 
there is only one arbitrary function, Ojy(r), and there is one free parameter less. The 
model is not only appealing because there are fewer degrees of freedom: only in a model 

| Note that while apparently this model is similar to that of Ref. [2D], it differs in the details. They fix 
their gauge A(r,t) — Ao(r) — r at the moment of recombination, while we fix it at the current time, 
and they define their free functions in terms of F(r) and k(r), while we use H(r) and £l(r). 
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with a homogeneous time to Big Bang the void corresponds to a perturbation in a FRW 
model at early times with no decaying mode [2T]. Ofcourse, if the void is not a very 
rare fluctuation in the density field, but created through a causal process in the early 
universe, this may not neccesarily be a constraint. 



2.5. Apparent acceleration in the light cone 

It is yet not clear, if what we observe as a function of redshift, in the form of luminosity 
distances to standard candles (e.g. Supernovae Type la), angular diameter distances to 
standard rulers (e.g. Baryon Acoustic Oscillations), the galaxy power spectrum, galaxy 
cluster counts, or other measures of the geometry and mass distribution of the universe, 
are due to modifications of gravity, an extra energy component, a cosmological constant, 
or simply a wrong interpretation of the underlying cosmological model. But as a whole 
they represent different possible explanations for the "Dark Energy" problem. One of 
the main observables, that will help decide between the different scenarios in the future, 
is the Hubble parameter H(z). Under the assumption that the correct background is a 
flat FRW cosmology we can write it [22] 



Hi j (z) 

' = (i + z) 3 n in + (i-n in )ex P 



j' +Z dlog(l + z')(l + w^(z')) 



Tj2 \ / 111 v --my — jr i or 1 /\ 1 -~ett \ 

- .,. - (49) 

-"in L Jl J 

where H in and f2 in are the expansion rate and matter density as observed at z = d§|. By 
taking the derivative we can write it as 



tl l^ log 

W eS ( Z > = -1 + 



%^-(l + 2) 3 fi i: 



3 dlog[l + *] 



(50) 



( 2 ) ■ 

where we have assumed that —^k (1 + z) z Vt- m > 0. The beauty of Eq. (1301) is that 

in 

if the observational data indeed is a manifestation of extra energy components, w^ L (z) 
has the usual interpretation of a dark energy equation of state, while in the case of 
modified gravity models, or the LTB model that we are considering in this paper, it 
can be interpreted as an empirical observational signature. In those case w^ L (z) is a 
function that captures the difference between the expansion rate that we measure, and 
the expansion rate that we ascribe to the observed matter density Q m (see also Eq. [37| 
for a correct definition of w in the case of an LTB universe). 

It is worthwhile pointing out, that even if there is not an accelerated expansion in 
the LTB models we are considering, because data is observed in the light cone, the total 
time derivative is [19] 



D__d_ cy/1 - grj d ^ d d_ 

Dt~ di A'{r,t) dr~ ~ dt ~ C ~dr~ ' ^ ' 

and an observer can measure an apparent acceleration when the light cone traverses 
the central inhomogeneity due to spatial gradients in either the matter density or the 

§ Alternatively, in [IS] Eq. (|50[> is written with derivatives of the scale factor a instead of the redshift 
z. We advocate using z, as it is an observable, in contrast to a(z). Notice that Hx(0,t) = i?i(0,f). 
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Figure 2. The apparent acceleration w^ L (z) inferred by applying Eq. (f5TT|) and FRW 
cosmology as the underlying model to describe the change in the Hubble parameter 
(either Ht(z) or Hl(z)) in the two best fit LTB models. 



expansion rate. In Fig. [5] we show w^ L (z) for two best fit models. Letting H = 
and H = Hl are both relevant, because different observations probe different expansion 
rates, i.e. as will be seen below, the Baryon Acoustic Oscillation signal depends partly 
on Hl, while supernova observations are only related to Ht through its dependence on 
di — (1 + z) 2 dA = (1 + z) 2 exp(J Hxdt). Interestingly, the variation and derivative of 
w^ L (z) is quite large in the best fit LTB models, showing that a precise low redshift 
supernova survey, such as the SDSS Supernova Survey [24J, sensitive to Ht, or a fine 
grained BAO survey such as PAU [2S], sensitive to H L , could rule out or reinforce 
the models in the near future. Conversely, if a disagreement between w as observed 
by Supernovae and w as observed through the BAOs is found, this could be a hint of 
inhomogeneous expansion rates. 

We can directly compute w^ L (z) in the limiting cases z = and z ^> 1 for 
asymptotically flat LTB models 

f "I + IjrSlW if * = 0, and H = H T 
= < -l + lj^k if - = 0, and H = H L , (52) 
[ if z > 1 

where we have used that the LTB metric converges asymptotically to a FRW metric 
giving dz = —da/ a. We see that to have w — 1/3 at low z implies either a significant 
negative gradient in H (r), or Q in ~ 1. 

2.6. Physical scales in the early universe 

Many bounds from observational cosmology, such as the sound horizon (a "standard 
ruler"), the CMB, and the big bang nucleosynthesis, are derived by considering scales 
and processes in the early universe, and are based on the implicit assumption of an 
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underlying FRW metric. To test LTB models against these observational data we have 
to connect distance scales, redshifts, and expansion rates in the early universe to those 
observed today. 

By construction, at high redshifts the LTB metric converges to a FRW metric, and 
the central void disappear (see Eq. flOl) . and physical results derived for FRW space-times 
still hold in the early universe, even though we are considering an LTB space-time. But 
starting from an approximately uniform universe at a high redshift z e in the LTB model, 
the expansion rate and matter density become gradually inhomogeneous, and a uniform 
comoving physical scale / in the early universe at z e , for example the sound horizon, is 
not uniform at some later redshift z. In particular the comoving size at t — to — ^bb (0) 
depends on how much relative expansion there has been at different positions since the 
formation of the uniform scale 

l(r(z))-l(r ) A (r^t(z e )) 

where to is the time now for the central observer, and is the radial coordinate of an 
observer very far away from the void. This is a consequence of defining the comoving 
physical scale as the scale measured at t . If instead we fixed the comoving length scales 
to be measured in the early universe at t(z e ), then indeed l{r(z)) would be independent 
of the observer position. The convenience of the above formula is, that the LTB models 
we consider are asymptotically FRW, and we can easily compute comoving scales at 
infinity. 

Normal relations determining early universe quantities, such as the redshift at 
matter-radiation equality, are essentially based on the Friedmann and conservation 
equations to relate cosmological parameters now to the parameters then 

H\z) = Hi [Q M (1 + zf + (1 - n M )(l + zf] (54) 
p(z)a 3 = p(0)a 3 , (55) 

where we have written it for a matter dominated universe. Eq. ( jTUl) is the LTB equivalent 
to Eq. (1541) . and since we are considering LTB metrics that asymptotically converge to 
FRW metrics the equation for light rays, or null geodesies, (115j) at high redshifts has the 
usual solution d(l + z) — —da/ a, and we can write an asymptotic version of Eq. ffTUl) 
that is valid at high redshifts 

H\z) = Hl s [n£(l + zf + (1 - n£)(l + zf] (56) 

where 

A (r(z e )) i 3/2 



H e g = Ho(z e 



(57) 



(1 + z e )A(r(z e ),t(z e ))_ 
and the asymptotic matter density is 

nf I = n M (r(z e )) = n out = i (58) 

In summary, any quantity at high redshifts in the LTB model can be computed with 
the usual formulas valid for a FRW metric, but using the matter density and Hubble 
constant given in Eqs. (l5"Tj)- (l58"l) . 
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3. Observational Data 

To assess the viability of the proposed models we have tested them against a set of 
current observational data. We divide the data into two classes: Constraining data sets, 
and prior data. The constraining data sets are actual measurements with errors, that 
are used to compute the likelihood of a given model, while the prior data, merely give 
ranges inside which the models should be. 

3.1. The Cosmic Microwave Background 

It is not a priori clear how to compute the spectrum of temperature anisotropies, without 
a full perturbation theory for the LTB models, but as shown in section 12.61 we can 
compute the comoving distance to the surface of large scattering, and the comoving size 
of the sound horizon at large distances. The ratio give the typical size, of the CMB 
temperature fluctuations on the sky, or equivalently the scale of the first peak, which is 
measured with exquisite precision by the WMAP satellite [26] 

CMB = ^7 rec { = 0.5952 ± 0.0021° = 0.010388 ± 0.00037. (59) 

The sound horizon r s (z Iec ) is calculated using the fitting formula provided by Eisenstein 
and Hu [21] with fl M — Q out = 1, h = H cS (see Eq. EZj). The physical baryon density 
in the early universe and the recombination redshift are fixed to their best fit WMAP3 
values QbH^ = 0.0223 and z Tec = 1089. The \ 2 from the CMB constraint is simply 

■-roc 1 1 

XCMB = 2 ( 60 ) 

6 

3.2. Baryon Acoustic Oscillations 

The BAO has been measured at different scales using a variety of techniques, and the 
feature has been detected in the 3D two-point correlation function [281 129] , the 3D power 
spectrum [3D], and the angular power spectrum [31J. In particular Percival et al [3DJ 
have combined the 2DF and SDSS large scale surveys to yield a measure of the BAO 
centred at two different redshifts, namely z = 0.2, and z = 0.35. 

In [3D] the power spectrum is calculated in a reference cosmology, and then the 
comoving distance scale is either dilated, using a fixed factor, or deformed using a 3 node 
spline fit. In principle one would need a full perturbation theory for LTB geometries to 
recalculate the 3D power spectrum, using our best fit model, but if we take into account 
that the observed galaxies are divided into redshift slices, and that the LTB universe 
at constant redshift behaves locally like a homogeneous FRW universe (e.g. Eq. (TTTJT) 
is a local analogy to the normal FRW Hubble rate equation), we can then, as a first 
approximation, relate the ACDM power spectrum calculated in [30J to the LTB power 
spectrum through a dilation. In Fig. [3] we show the fractional difference between a 
simple dilation of the comoving scale, and a full modelling around the relevant redshifts 
z = 0.2 and = 0.35. The dilation is an excellent approximation, the difference being less 
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Redshift Redshift 

Figure 3. The fractional difference in the comoving distance used in the reference 
model in [3D] compared to the best fit GBH model (red, full line) and constrained 
model (blue, dashed line). To the left (right) is shown the the dilation at z = 0.2 
(z = 0.35), and in parenthesis in the legend is indicated the dilation. 



than 2% over the relevant redshift range. We stress, however, that a comprehensive test 
of LTB models against large scale structure data has to await the development of the 
linear perturbation theory for LTB space-times, an approach that is outside the scope 
of this paper [32] . 

The relevant BAO quantity to test against is the observed peak (in the case of 
real space) or wiggle wavelength (in the case of Fourier space), compared to the size 
of the sound horizon in the reference model used in [30], which can be interpreted 
approximately as the projected size of the BAO on the sky for structures at a given 
redshift. The ratio is given as 

0bao(z) = $L. (61) 
D v {z) 

As discussed in Section 12.61 the comoving size of the sound horizon r s (z) is a 
function of redshift, due to the inhomogeneous nature of the expansion since the surface 
of last scattering A(r(z), t(z))/A(r(z Tec ), t(z Tec )). Because we detect the BAO in the 3D 
distribution of galaxies, and not in the 2D projection on the sky, Dy is not just the 
comoving distance, but rather a combination of longitudinal and transversal distances. 
In LTB space-times the longitudinal and transversal Hubble parameters, Hl and Ht, can 
differ significantly, see Fig. 1, and therefore it is important to use the correct longitudinal 
expansion rate: 

-,1/3 

(62) 



D v (z) 



d\(z)(l + zf ' 



Hl{z) 
Percival et al [30] find 

# BAO (0.2) = 0.1980 ± 0.0058 #bao(0.35) = 0.1094 ± 0.0033 (63) 

with a 39% correlation, and we use these two measurements as our BAO data set giving 
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the x 



with 



2 



Xbao = 0BAofc) - r s { Zi )/D v { Zi )] [9 BAO { Zj ) - r s ( Zj )/D v ( Zj )} (64) 



T-l 



35059 -24031 
-24031 108300 



(65) 



3.3. Type la Supernovae 



We use the Type la Supernovae compilation by Davis et al. [33], which is a compilation 
of 192 SNIa consisting of 45 SNe from a nearby SNIa sample [34J, 57 SNLS [35J and 
60 ESSENCE [36J intermediate redshift SNe and 30 high redshift "Gold" SNe [37], 
with internally consistent magnitude offsets. The supernovae span the redshift range 
z = 0.01 — 1.7, and we use the magnitude residuals \i to constrain the LTB models. 
The residual /x, and apparent and absolute magnitudes m and M are related to the 
luminosity distance 4 as 

d L 



H = m — M = 5 log 



10 



IMpc 



+ 25 (66) 



The exact absolute magnitude of a SNIa is unknown, and we include an arbitrary offset 
/io, when calculating \ 2 for the model fit to the observed SNIa 

V [/^-(/^(*) +/*>)] 2 ( , 7) 
~ erf 

where /i is determined by minimising XsNia- Z 1 nas a logarithmic dependence on 4, 
and the zero point /io is degenerate with the local overall scale of the expansion rate 
(H in in the GBH model, and H in the constrained model, see Figs. [HlandE]). 

Even though the Supernova data set is by far the largest of the three, the error bars 
on individual SNe are large, and internally there is a large scatter, as can be appreciated 
in Fig. HI where the residuals are compared to the best fit GBH LTB models, and the 
open CDM and ACDM FRW models. Note that the predicted curves for the best fit 
GBH model and the best fit ACDM model start to deviate significantly beyond redshift 
z = 1, and therefore it would be extremely useful to have a complete Supernovae data 
set at high redshifts, which could help discard one of the models against the other one. 

The three data sets discussed above are obtained through different means, and relate 
to different physical observations, and hence they are not correlated in any way. This 
independence allow us to find the total likelihood by simply multiplying the individual 
likelihood functions. 



3-4- Priors 

Even though we do not include them in the likelihood analysis, we still require our 
models to obey three additional priors: Two of them concern the local universe, and as 
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Figure 4. Apparent magnitude residuals for the two best fit LTB models, standard 
open CDM, and the best fit ACDM FRW models compared to Type la Supernovae 
data. 



such are priors at z — 0. They are the observed lower age limit on globular clusters in 
the Milky Way of 11.2 Gyr |38j (2-cr limit) and the HST key project [35] measure 
of the local value for the Hubble parameter H m = 72 ± 8 (1-er limit). The third 
prior is the gas fraction as observed in clusters of galaxies. This is a very powerful 
observation for limiting alternative models: Clusters of galaxies sit at the bottom of 
deep gravitational potential wells, and supposedly the gas fraction is representative 
for the universe as a whole, because neither gas nor dark matter can escape out 
of the potential well. This universal prior can be compared to the gas fraction we 
deduce far away from the void, by combining the WMAP satellite and Big Bang 
Nucleosynthesis bound on the physical baryon density uj^ = QbH% s = 0.0223 with the 
physical matter density at infinity uj m = Q out H^ s = H^ s . The current observational 
limits are / gas = uJb/oj m = 0.1104 ± 0.0016 ± 0.1 (random+systematic) [10], while for 
our best fit models (see Table ED we find / gas = 0.127 — 0.134 in agreement at 2-cr with 
observations. 



Model 


Ho H m H out 




ro Ar 


units 


100 km s" 1 Mpc" 1 




Gpc r 


GBH 


0.5 - 0.85 0.30- 0.70 


0.05-0.35 


0.3-4.5 0.1-0.9 


Constrained 


0.50 - 0.95 0.4 - 0.89 0.33 - 0.63 


0.05-0.35 


0.5-4.5 0.1-0.9 



Table 1. Priors used when scanning the parameters of the two models. In the 
constrained model Hq is only a pre factor for Hq(t) and the span of H m and flout 
are derived from the priors on VL m and H - 




Figure 5. To the left/right the CMB spectrum for the best fit GBH/constrained-GBH 
model (red line) compared to WMAP3 data (diamonds) 42J, and including cosmic 
variance (grey shading). 



4. Analysis and results 

To test the full and the constrained GBH LTB model, we have performed a parameter 
scan over the models and for each set calculated the \ 2 ■ The priors are given in 
table [D, and have been chosen to encompass the best fit 2-a limits, except where large 
degeneracies exist, and also to be reasonable, taking into account the HST key project 
[39] . and acceptable matter densities. 

4-1. The best fit models 

In table [5] are given the best fit for the two models. It is interesting to notice that both 
have very similar values, and in fact this seems to hint that current data prefer the 
simpler constrained model with a homogeneous big bang. Both best fit models have 
local Hubble rates on the low side but still in agreement with the HST project, and the 
local time to Big Bang is well inside the limits given by globular clusters. 

The best fit models both give an excellent luminosity redshift distance relation, that 
are in as good an accordance with current Type la SNe as the ACDM model (see Fig. 
H]). This comes as no surprise: LTB models can be constructed that fit any luminosity- 
redshift relation [H]. Our model fit is done though under the simultaneous constraints 
of the other probes, and as such is more constrained. In the transition zone between the 
void and the surrounding Einstein-de Sitter space (at r ~ r ) there is a significant (up 
to ~10%) difference between Ht{t), which is related to d^, important for supernova 
observations, and Hl{v), important for the longitudinal part of the Baryon Acoustic 
Oscillations (see Fig. []]). This difference marks a fundamental observational signature 
between LTB and FRW models. A set of very well observed Type la SNe at z 0.1 — 0.2, 
such as the SDSS II SN survey [21], together with the BAO observations already done 
at similar redshifts, will put strain on either model. 

Because we currently do not have a full perturbation theory for LTB space- 
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times, we were not able to make a full likelihood analysis comparing our model to 
all the WMAP data. Nonetheless, to get an idea of how bad the models fit the full 
body of WMAP temperature anisotropy observations, we have calculated a standard 
temperature anisotropy spectrum using as input H eS and Q = Q out = 1 (see Fig. 
Even though we fix the physical baryon density, and only fit the first peak of the 
CMB, the obtained model is not too bad, and it is reassuring that other people have 
proposed Einstein-de Sitter models that do fit the WMAP3 CMB observations using 
a non standard primordial spectrum and a hot neutrino component (e.g. |43j). Very 
recently an LTB model similar to ours was proposed jH] and the authors managed to 
make a reasonable fit to the full WMAP TT and TE data by adding a running of the 
tilt in the primordial spectrum. In a more complete analysis at low Z-values and large 
angles one could expect an effect similar to that of a cosmo logical constant, since we 
have a non-trivial curvature at low redshifts, and hence an even better fit that do not 
have to rely on ad hoc features or a large running of the tilt in the primordial power 
spectrum may indeed be possible. 

4-2. Likelihood contours and degeneracies 

We have marginalised over different dimensions in parameter space by integrating over 
the likelihood, given as £ oc exp(— % 2 /2). The marginalised l-o~ and 2-cr likelihoods for 
the individual data sets, and also the 3-cr limits for the combined data sets are shown 
in Figs. M and [6j 

In the normal ACDM model, if cj& is fixed, there is a well known strong degeneracy 
between Q m and i^ojlH] for a given size of the sound horizon on the sky 9 a- In our model 
the relevant Q m for the CMB is Q out = 1, that fixes the size of the sound horizon, and 
we have also fixed u>b- Then the size of 9a depends essentially on H or H in , which is 
reflected in the likelihood constraint from CMB on Hubble rates as seen in the figures. 
Nonetheless we can to some extent change dA by introducing curvature, allowing us to 
choose a higher value for or Hq, either by having a large void size, ro, or making 
the void very underdense Q[ n <C 1. This can be seen by the widening of the 2-cr limit in 
the H-vq plot for large values of ro, and the asymmetric 2-cr errors in the H-Qi n plots. 

f2; n and ro are the major parameters determining the luminosity-redshift relation, 
and are hence constrained by the Type la SNe and the BAO. While only the relative 
value of Hubble rate play a role for the SNe, the BAO does limit H, and there is some 
strain between the BAO and the Type la SNe, as seen in the fi in -r plots. 

An obvious degeneracy is that of ro and Ar/ro because the effect of a larger void 
can also be obtained by making a smoother, and hence broader, transition. Current 
data does not have any sensitivity to H ont , because it mainly affects the Hubble rate at 
very high redshift, where no good observational data exist. The relative transition width 
Ar/r is also not very well constrained. The only thing we can deduce, in agreement 
with the good fit to the data given by the ACDM model, is that no sudden transition 
is allowed. The lower limit on Ar/r is mainly limited by the Type la SN data, and to 
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Figure 6. Likelihoods for the GBH model: The likelihood for the combined data 
set is shown in yellow with 1-, 2-, and 3-er contours, while the individual SNIa, BAO, 
and CMB data sets are shown in blue, purple, and green respectively with 1- and 2-a 
contours. 

a lesser extent the BAO. 

It should be stressed though, that even though our LTB models give very good fits 
with x 2 that are comparable to that of the ACDM model, current data do put significant 
constraints on the models, and they will probably be challenged by new observational 
data in the near future, and can be falsified. 



Model 


Hq H m H out H cS 


f2; n 


r Ar 


^BB 


units 


100 km s- 1 Mpc- 1 




Gpc r 


Gyr 


GBH 


0.58±0.03 0.49±0.2 0.43 


0.13±0.06 


2.3±0.9 0.62(>0.20) 


14.8 


Constrained 


0.64±0.03 0.56 0.43 0.42 


0.13±0.06 


2.5±0.7 0.64(>0.21) 


15.3 



Table 2. Best fit values with 2-a error bars for the two models. The likelihood 
contours are not closed for Ar, and only a lower limit can be given. In the GBH model 
H out is unconstrained. For the central values of the other four parameters H out = 0.49 
minimises \ 2 ■ Notice that naturally the best fit GBH model and the constrained model 
give similar best fit values, and error bars. Ie. among all the different GBH models a 
model with a homogeneous Big Bang is preferred. 
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Figure 7. . . . Figure [5] continued. 
4-3. Bayesian analysis 

In this section we would like to find out whether a homogeneous FRW model of the 
universe (including the accelerated expansion in terms of a vacuum energy with constant 
equation of state w) can be used with confidence when analysing present cosmological 
data, or should we rather be more general and assume an inhomogeneous LTB model 
of the universe? 

The standard frequentist analysis of parameter estimation, given a set of data, is 
not very useful for model selection, since it is difficult to compare models with different 
number of parameters. For a discussion about probability theory and model selection 
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Figure 8. Likelihoods for the constrained model: The likelihood for the combined 
data set is shown in yellow with 1-, 2-, and 3-er contours, while the individual SNIa, 
BAO, and CMB data sets are shown in blue, purple, and green respectively with 1- 
and 2-cr contours. 



see Refs. [I6j H7J HH H9]. For instance, the usual method of comparing minimum \ 2 per 
effective degree of freedom normally misses the point and is not very decisive. Other 
methods to decide which model gives the best description, given the data, include various 
Information Criteria, e.g. Akaike joO] and Bayesian [ST], which use more or less ad hoc 
formulae without much justification and normally do not compare well among eachother. 
However, in the last few years there has been a flourishing of several independent 
analysis based on the Bayesian evidence associated with a given likelihood and a given 
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cosmological model, within some given priors, both theoretical and observational, see 

e.g. Refs. [52H53H5HE5]. 

The Bayesian evidence is based on Bayes theorem, which relates the posterior 
distribution V(8, JA\D) for the parameters 9 of the model M. given the data D, in 
terms of the likelihood distribution function C(D\6, M) within a given set of priors 
n(6,M), 

where E is the Bayesian evidence, that is the average likelihood over the priors, 

E(D\M)= I d6 £(D\6,M)ir{6,M). (69) 



The prior distribution functions contain all the information about the parameters 
before observing the data, e.g. our theoretical prejudices, our physical intuition about 
the model, the ranges of parameters obtained from previous experiments, etc. The 
Bayesian evidence is very useful because it allows a comparison of models within a 
complete set M.i=\...N (in our case, N=2). We can compute the posterior probability 
for each hypothesis (model) given the data D using again Bayes theorem, V(M.i\D) cx 
E(D\ M.i) ir(J\4i), where E(D\M.i) is the evidence of the data under model Mli and 
7r(Aii) is the prior probability of the i-th model before we see the data, usually taken 
to be identical, i.e. 7r(A^j=i...jv) = l/N. Finally, the ratio of the evidences for two 
competing models is called the Bayes factor, 

_ £(D|M) , . 

Bij = EWMj) ■ (70) 
This expression provides a mathematical representation of Occam's razor, because more 
complex models tend to be less predictive, lowering their average likelihood (within the 
priors) in comparison with simpler, more predictive models. Complex models can only 
be favoured if they are able to provide a significantly improved fit to the data. In simple 
cases where different models give vastly different maximum likelihoods there is no need 
to employ model selection techniques because they provide only minor corrections to 
the standard inference, but they are essential when the difference between maximum 
likelihoods is only marginal, as will be the case at hand. The Bayes factor ( J7D1 is then 
used to give evidence of (i.e. favour) the model M.i against the model M.j using the so 
called Jeffreys' scale, a particular interpretation of the Bayes factor which strengthens its 
veridic roughly each time the logarithm lni?^ increases by one unit, from (undecisive) 
to greater than 5 (strongly ruled out). 

Unfortunately, the computation of the Bayesian evidence (1691) is rather involved 
and typically requires extensive computational power, unless the number of parameters 
is significantly reduced. When the likelihood of the data given the model parameters 
is a single isolated peak, far from the edges of the prior ranges, then there is a simple 
approximation to the logarithm of the Bayesian evidence, 

n 

ln£ = ln£ max -lnA-^lnA0i, (71) 
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where A is the normalisation of the likelihood, and A6i = bi—ai is the range of parameter 
6i G [dj, bi], i — 1 . . . n. Moreover, for the case of a Gaussian likelihood, 



C(9) — A exp 



2 



we find A = (27r) -n / 2 /Vdet C, where C is the covariance matrix and Xi — 9i — 6~i. It is 
clear that whenever the prior ranges are too big for the likelihood, the Bayesian evidence 
is penalised. Moreover, the more parameters there are, the larger the difference in InE, 
and the larger the logarithm of the Bayes factor, as expected. 

We will now apply this discussion to the case at hand. We will compare the standard 
flat ACDM model with a constant equation of state parameter, w, while fixing all other 
cosmo logical parameters but Q M - Thus our FRW model has here 2 parameters, which 
we believe are the most constraining, since others like the local rate of expansion do not 
provide significant extra information. We performed a calculation of the minimum x 2 
within a grid of 5000 models, where w G [—2.0, 0.0] and Vl M G [0.0, 0.6]. The maximum 
likelihood corresponds to a x^ACDM) = 197.05 for w = -1.005 and Vt M = 0.276. 

On the other hand, the constrained GBH model has 4 independent parameters, 
the local rate of expansion Hq = 100 h km/s/Mpc, the local matter fraction Q m = 
^m(0), the transition distance r , and the width Ar/r , providing a grid of several 
million models within the ranges h G [0.50,0.95], Q m G [0.05,0.35], r G [0.5,4.5] 
and Ar/r G [0.1,0.9]. We found that the maximum likelihood corresponds to a 
X 2 min (GBH) = 197.845 for h = 0.659, tt m = 0.124, r = 2.47 Gpc and Ar/r = 0.638. 

At face value it seems that the inhomogeneous model provides as good a fit to the 
data as the FRW one. If we compute the usual minimum x 2 P er effective number of 
degrees of freedom (i.e. the number of data points minus the number of parameters), 
we find 

X 2 nin (ACDM)/d.o.f. = 1.021 , xl m (GBH)/d.o.f. = 1.036 , (72) 

so that both models seem excellent descriptions of the data, the first one being slightly 
better. However, there are other indicators more appropriate for model comparison, like 
the (corrected) Akaike Information Criterion (AIC), computed as 

AIC = ;d m + 2k + ^'-l ' ^ 

where k is the number of parameters and is the number of data points. In our case this 
gives AIC(ACDM) = 201.1, while AIC(GBH) = 206.0, which would clearly favour 
the homogeneous FRW model. On the other hand, if we choose to compare models with 
the Bayesian Information Criterion (BIC), computed as 

BIC = xLn + klnN, (74) 

we find BIC(ACDM) = 207.6, while BIC{GBH) = 218.9, which would very strongly 
favour the homogeneous FRW model. Clearly, neither method gives a good assessment 
for choosing among models. This is the reason why the Bayesian evidence has been used 
recently in the context of model comparison. 
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If we compute the Bayes factor (170j) by performing the integral of the likelihood 
over the priors, $M) , we find 

In E(ACDM) = -103.1, In E(GBH) = -106.7, (75) 

and therefore the logarithm of the Bayes factor is ln5 12 = 3.6, which clearly favours 
the homogeneous FRW model against the GBH-LTB model. It seems that the bayesian 
evidence method discards significantly, but not very strongly, the inhomogeneous model 
versus the usual FRW model. It is possible that, in the future, better data sets 
and stronger priors on cosmological parameters may discard once and for all the 
inhomogeneous model. For the moment, a local void with a size of several Gpc, with 
matter well below average, and a local rate of expansion of 71 km/s/Mpc, can account 
for both the distant (CMB), intermediate (SNIa) and local (BAO) data sets. 

5. Discussion and Conclusions 

We have shown that present observations do not exclude the possibility that we live close 
to the center of a large void. This is a appealing possibility which effectively gets rid of 
the necessity to introduce an ad hoc cosmological constant in our model of the universe. 
Moreover, it is consistent with early universe cosmology in terms of inflationary initial 
conditions for the origin of large scale structures. Perhaps these voids arise due to large 
non-perturbative inhomogeneities associated with the stochastic nature of the inflaton 
evolution [15], or due to large non-gaussianities in the primordial spectrum coming from 
inflation that could arise due to phase transitions or in multifield inflation [56, 57J. In 
fact, we already know there can be other voids in our local patch of the universe, with 
a Gpc scale, as exemplified by the observed cold spot in the CMB [T4] . 

We have analysed the likelihood of such an interpretation of the present acceleration 
of the universe, using data from the Cosmic Microwave Background (on large scales), 
Supernova la (at intermediate scales), Baryon Acoustic Oscillations, the present age 
and the local rate of expansion (at small scales). All the data seems to be consistent at 
the 95% confidence level with a local void of size around 2.5 Gpc, within an Einstein-de 
Sitter universe on large scales, without the need to introduce a cosmological constant. 
The apparent acceleration of the universe can be interpreted here as due to the curved 
path of photons in this locally open universe. We should however take into account that 
a possible displacement from the centre of the void will induce a dipole in the CMB 
[SHI EH] • If this (geometrically induced) dipole is not cancelled by a peculiar velocity 
towards the centre, to be in agreement with the opbserved CMB dipole, one can show 
[59] that we cannot be displaced more than approximately 100 Mpc from the centre. 
From a probabilistic point of view, this is a generic argument against Gpc sized voids. 

We have performed a Bayesian analysis in order to compare two competing models, 
the predominant ACDM model and our GBH inhomogeneous model. While the usual 
frequentists analysis does not discard the GBH model against the ACDM model, there 
seems to be strong but not decisive (bayesian) evidence against the GBH model. It 
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is possible that in the near future, with much better cosmological data on large scale 
structures and an extended set of supernovae at intermediate and high redshift, we 
may be able to constrain and definitely rule out the inhomogeneous LTB model. We 
should also mention that the data seems to favour a homogeneous Big Bang since the 
constrained GBH model gives the same likelihood contours and minimum \ 2 than the 
unconstrained model, while having one parameter less. 

At the moment we are studying the effect that a generic LTB model has on the 
growth of structure in order to constrain further the GBH model with data from ISW- 
LSS correlations and the Lyman-a forest within the Alcock-Paczynski-test analysis [32J. 

In conclusion, we cannot discard that we live in an inhomogeneous local void 
within an asymptotically Einstein-de Sitter universe. The possibility that we have 
misinterpreted the present acceleration and that the cosmological constant is nothing but 
a mirage has been addressed recently [60, EH E2J [63] . We have added to the discussion 
the comparison with a large, albeit incomplete, set of cosmological observations, and 
the bayesian analysis appropriate for model selection. We hope in the future to provide 
further constraints on the model and possibly rule it out. 

Acknowledgements 

We thank Andrei Linde, Andy Albrecht, David Valls-Gabaud, Sysky Rasanen, and Alex 
Kusenko for generous comments and suggestions. We also thank the Danish Centre of 
Scientific Computing (DCSC) for granting the computer resources used. JGB would 
like to thank the Kavli Institute for Theoretical Physics for hospitality during the last 
stages of the work, supported in part by the National Science Foundation under Grant 
No. PHY05-51164. We also acknowledge financial support from the Spanish Research 
Ministry (M.E.C.), under the contract FPA2006-05807. 

References 

[1] http : // www . darkenergy survey .org. 

[2] http : //www. aip . org/web_bin/pt/vol-58/ iss-8/p53 . shtml 

[3] George Lemaitre. " The Expanding Universe" . Gen. Rel. Grav., 29:641-680, 1997. 
[4] Richard C. Tolman. "Effect of imhomogencity on cosmological models". Proc. Nat. Acad. Sci., 
20:169-176, 1934. 

[5] H. Bondi. " Spherically symmetrical models in general relativity" . Mon. Not. Roy. Astron. Soc, 
107:410-425, 1947. 

[6] Idit Zehavi, Adam G. Riess, Robert P. Kirshner, and Avishai Dekel. "A Local Hubble Bubble 

from SNe la?". Astrophys. J., 503:483, 1998, astro-ph/9802252. 
[7] Kenji Tomita. "Anisotropy of the Hubble Constant in a Cosmological Model with a Local Void 

on Scales of 200 Mpc". 2000, astro-ph/0005031. 
[8] Kenji Tomita. "A Local Void and the Accelerating Universe". Mon. Not. Roy. Astron. Soc, 

326:287, 2001, astro-ph/0011484. 
[9] Kenji Tomita. "Analyses of Type la Supernova Data in Cosmological Models with a Local Void". 

Prog. Theor. Phys., 106:929-939, 2001, astro-ph/0104141. 



Confronting Lemaitre-Tolman-Bondi models with Observational Cosmology 



26 



William J. Frith, G. S. Busswell, R. Fong, N. Metcalfe, and T. Shanks. "The Local Hole in the 
Galaxy Distribution: Evidence from 2MASS". Mon. Not. Roy. Astron. Soc, 345:1049, 2003, 
astro-ph/0302331. 

G. S. Busswell et al. "The Local Hole in the Galaxy Distribution: New Optical Evidence". Mon. 

Not. Roy. Astron. Soc, 354:991, 2004, astro-ph/0302330. 
Patricio Vielva, E. Martinez-Gonzalez, R. B. Barreiro, J. L. Sanz, and L. Cayon. "Detection of 

non-Gaussianity in the WMAP 1-year data using spherical wavelets". Astrophys. J., 609:22-34, 

2004, astro-ph/0310273. 

M. Cruz, M. Tucci, E. Martinez-Gonzalez, and P. Vielva. "The non-Gaussian Cold Spot in WMAP: 
significance, morphology and foreground contribution" . Mon. Not. Roy. Astron. Soc, 369:57-67, 
2006, astro-ph/0601427. 

Marcos Cruz, L. Cayon, E. Martinez- Gonzalez, P. Vielva, and J. Jin. "The non-Gaussian Cold 

Spot in the 3-year WMAP data". Astrophys. J., 655:11-20, 2007, astro-ph/0603859. 
Andrei D. Linde, Dmitri A. Linde, and Arthur Mezhlumian. "Do we live in the center of the 

world?". Phys. Lett., B345:203-210, 1995, hcp-th/9411111. 
Andrei D. Linde, Dmitri A. Linde, and Arthur Mezhlumian. "From the Big Bang theory to the 

theory of a stationary universe". Phys. Rev., D49:1783-1826, 1994, gr-qc/9306035. 
Juan Garcia-Bellido, Andrei D. Linde, and Dmitri A. Linde. "Fluctuations of the gravitational 

constant in the inflationary Brans-Dicke cosmology". Phys. Rev., D50:730-750, 1994, astro- 

ph/9312039. 

Kari Enqvist and Teppo Mattsson. "The effect of inhomogeneous expansion on the supernova 

observations". JCAP, 0702:019, 2007, astro-ph/0609120. 
Kari Enqvist. "Lemaitre-Tolman-Bondi model and accelerating expansion". 2007, 

arXiv:0709.2044 [astro-ph]. 
Havard Alnes, Morad Amarzguioui, and Oyvind Gron. "An inhomogeneous alternative to dark 

energy?". Phys. Rev., D73:083519, 2006, astro-ph/0512006. 
Charles Hellaby and Andrzcj Krasinski. Alternative methods of describing structure formation in 

the Lemaitre-Tolman model. Phys. Rev., D73:023518, 2006, gr-qc/0510093. 
Eric V. Linder and Adrian Jenkins. "Cosmic Structure and Dark Energy". Mon. Not. Roy. 

Astron. Soc., 346:573, 2003, astro-ph/0305286. 
Eric V. Linder. "Cosmic growth history and expansion history". Phys. Rev., D72:043529, 2005, 

astro-ph/0507263. 



http : //sdssdp47 . f nal . gov/sdsssn/sdsssn . html 



http : //www . ice . csic . es/pau. 

D. N. Spergel et al. "Wilkinson Microwave Anisotropy Probe (WMAP) three year results: 

Implications for cosmology". Astrophys. J. Suppl, 170:377, 2007, astro-ph/0603449. 
Daniel J. Eisenstein and Wayne Hu. "Baryonic Features in the Matter Transfer Function". 

Astrophys. J., 496:605, 1998, astro-ph/9709112. 
Daniel J. Eisenstein et al. "Detection of the Baryon Acoustic Peak in the Large-Scale Correlation 

Function of SDSS Luminous Red Galaxies" . Astrophys. J., 633:560-574, 2005, astro-ph/0501171. 
Teppei Okumura et al. "Large-Scale Anisotropic Correlation Function of SDSS Luminous Red 

Galaxies". 2007, arXiv:0711.3640 [astro-ph]. 
Will J. Percival et al. "Measuring the Baryon Acoustic Oscillation scale using the SDSS and 

2dFGRS". 2007, arXiv:0705.3323 [astro-ph]. 
Nikhil Padmanabhan et al. " The Clustering of Luminous Red Galaxies in the Sloan Digital Sky 

Survey Imaging Data". Mon. Not. Roy. Astron. Soc., 378:852-872, 2007, astro-ph/0605302. 
Alicia Bueno, Juan Garcia-Bellido, and Troels Haugboelle. "Linear perturbation theory in LTB 

models" . in preparation. 
Tamara M. Davis et al. "Scrutinizing exotic cosmological models using ESSENCE supernova data 

combined with other cosmological probes". Astrophys. J., 666:716, 2007, astro-ph/ 0701510. 
Saurabh Jha, Adam G. Riess, and Robert P. Kirshner. "Improved Distances to Type la Supernovae 



Confronting Lemaitre-Tolman-Bondi models with Observational Cosmology 



27 



with Multicolor Light Curve Shapes: MLCS2k2". Astrophys. J., 659:122-148, 2007, astro- 
ph/0612666. 

Pierre Astier et al. "The Supernova Legacy Survey: Measurement of Om, and w from the 
First Year Data Set". Astron. Astrophys., 447:31-48, 2006, astro-ph/0510447. 

W. Michael Wood-Vasey et al. "Observational Constraints on the Nature of the Dark Energy: 
First Cosmological Results from the ESSENCE Supernova Survey". Astrophys. J., 666:694, 
2007, astro-ph/0701041. 

Adam G. Riess et al. "New Hubble Space Telescope Discoveries of Type la Supernovae at z > 1: 
Narrowing Constraints on the Early Behavior of Dark Energy". Astrophys. J., 659:98, 2006, 
astro-ph/0611572. 

Lawrence M. Krauss and Brian Chaboyer. "Age Estimates of Globular Clusters in the Milky Way: 

Constraints on Cosmology". Science, 299:65-70, 2003. 
W. L. Freedman et al. "Final Results from the Hubble Space Telescope Key Project to Measure 

the Hubble Constant". Astrophys. J., 553:47-72, 2001, astro-ph/0012376. 
S. W. Allen et al. "Improved constraints on dark energy from Chandra X-ray observations of the 

largest relaxed galaxy clusters". 2007, arXiv:0706.0033 [astro-ph]. 
Marie-Noelle Celerier. "Do we really see a cosmological constant in the supernovae data?" . Astron. 

Astrophys., 353:63-71, 2000, astro-ph/9907206. 

G. Hinshaw et al. "Three- year Wilkinson Microwave Anisotropy Probe (WMAP) observations: 
Temperature analysis". Astrophys. J. SuppL, 170:288, 2007, astro-ph/0603451. 

Paul Hunt and Subir Sarkar. "Multiple inflation and the WMAP 'glitches' II. Data analysis and 
cosmological parameter extraction". 2007, arXiv: 0706. 2443 [astro-ph]. 

Stcphon Alexander, Tirthabir Biswas, Alessio Notari, and Deepak Vaid. Local Void vs Dark 
Energy: Confrontation with WMAP and Type la Supernovae. 2007, arXiv:0712.0370 [astro- 
ph]. 

L. Page et al. First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Interpretation of the TT and TE Angular Power Spectrum Peaks. Astrophys. J. SuppL, 148:233, 
2003, astro-ph/0302220. 

H. Jeffreys. " The theory of probability" . Oxford U. P. (1998). 

E.T. Jaynes. "Probability Theory: the Logic of Science" . Cambridge U.P. (2003). 
D.J.C. Mackay. "Information theory, inference and learning algorithms" . Cambridge U.P. (2003). 
G. D'Agostini. "Bayesian reasoning in data analysis: A critical introduction". World Scientific 
(2003). 

Hirotugu Akaike. "A new look at the statistical model identification". IEEE Transactions on 

Automatic Control, 16 (6):716-723, 1974. 
G. Schwarz. "Estimating the dimension of a model". Annals of Statistics, 6 (2):461-464, 1978. 
Andrew R. Liddlc. " How many cosmological parameters?" . Mon. Not. Roy. Astron. Soc, 35LL49- 

L53, 2004, 

Maria Beltran, Juan Garcia-Bcllido, Julicn Lcsgourgucs, Andrew R Liddle, and Anze Slosar. 
"Bayesian model selection and isocurvature perturbations". Phys. Rev., D7L063532, 2005, 
astro-ph/0501477. 

David Parkinson, Pia Mukherjee, and Andrew R Liddlc. "A Bayesian model selection analysis of 

WMAP3". Phys. Rev., D73T23523, 2006, astro-ph/0605003. 
Roberto Trotta. " Applications of Bayesian model selection to cosmological parameters" . Mon. 

Not. Roy. Astron. Soc, 378:72-82, 2007, astro-ph/0504022. 
L. A. Kofman and Andrei D. Linde. Generation of Density Perturbations in the Inflationary 

Cosmology. Nucl. Phys., B282:555, 1987. 
Xingang Chen, Richard Easther, and Eugene A. Lim. Generation and Characterization of Large 

Non-Gaussianities in Single Field Inflation. 2008, arXiv:0801.3295 [astro-ph]. 
T. Paczynski, B. and Piran. A dipole moment of the microwave background as a cosmological 

effect. Astrophys. J., 364:341-348, December 1990. 



Confronting Lemaitre-Tolman-Bondi models with Observational Cosmology 



28 



[59] Havard Alnes and Morad Amarzguioui. "CMB anisotropies seen by an off-center observer 
in a spherically symmetric inhomogeneous universe". Phys. Rev., D74:103520, 2006, astro- 
ph/0607334. 

[60] Syksy Rasanen. Dark energy from backreaction. JCAP, 0402:003, 2004, astro-ph/0311257. 
[61] Teppo Mattsson. " Dark energy as a mirage" . 2007, arXiv:071 1.4264 [astro-ph] . 
[62] David L. Wiltshire. Dark energy without dark energy. 2007, arXiv:0712. 3984 [astro-ph]. 
[63] Nan Li, Marina Seikel, and Dominik J. Schwarz. "Is dark energy an effect of averaging?". 2008, 
arXiv:0801.3420 [astro-ph]. 



